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0^ ' Absract. Minimum total energy calculations, which account for both electron-lattice and electron- 

electron interactions in conjugated polymers are performed for chains with up to eight carbon atoms. These 
O ' calculations are motivated in part by recent experimental results on the spectroscopy of polyenes and con- 

jugated polymers and shed light on the longstanding question of the relative importance of electron-lattice 
vs. electron-electron interactions in determining the properties of these systems. 
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A large amount of experimental evidence [1] regarding conjugated polymers can be understood in terms 
of independent electron theories that account for electron-lattice (e-i) coupling and cr-bond compressibil- 
ity [2-3]. However, there exists a considerable body of spectroscopic results, concerning especially the order- 
f'S i ing of excited states [4-9], which cannot be explained without invoking electron-electron (e-e) correlations. 

CNJ ■ Since such different experimental results are usually rationalized in terms of models which describe ade- 

quately only either the e-J or the e-e interaction, different groups of researchers have been led to emphasize 
• in these systems the importance of one of the two effects at the expense of the other. 



In this letter the results of a set of minimum total energy calculations which fully include both interac- 
tions are presented. There are several reasons to pursue this goal. From a theoretical standpoint it is natural 
to assume that the transfer integrals depend on the distance between carbon sites and that there is an energy 
cost involved in stretching a carbon-carbon bond [2] . It is also not surprising to find manifestations of e-e 
pH ■ interactions which are not accounted for by models implying complete screening such as those of ref. [2-3]. 

Q , On the experimental side there is a growing amount of evidence indicating that the ordering of excited states 

O ' depends on the specific polymer and in some instances it appears that different probing techniques lead to 

different results in this regard [4-9]. In particular, recent observations in short thiophene oligomers [6] and 
in poly(p-phenylene-vinylene) [7] show that in these systems the ordering of the two lowest excited states is 
reversed compared to that observed in polyenes [4-5] . Within the context of the SSH model [3] it is natural 
' to interpret this reversal in terms of lack of ground state degeneracy in the systems of ref. [6-7] . This is be- 

cause, upon excitation of one electron from the highest occupied to the lowest unoccupied molecular orbital, 
lack of ground state degeneracy leads to two separate bipolaron levels as opposed to a pair of degenerate 
soliton levels. E-e repulsion favors the 2 ^Ag over the 1 level [10], and the results of ref. [6-7] suggest 
that this latter effect is not strong enough to overcome the energy difference between bipolaron levels in 
these systems. These qualitative considerations hint to the possibility that important physical effects may 
be overlooked if the spectroscopic results are interpreted without fully accounting for e~l interactions. 

Minimum total energy calculations based on the SSH description were presented in ref. [10-11]. For 
a iraras-polyacetylene chain with N carbon atoms, the starting point is the hamiltonian for the vr electron 
system: 

^ = - X! (to~a{un+l~Un)j ■ {cl+i^sCn,s + 4 ,,C„+i_^) + ^ ^(u„+i-U„)^ (1) 

Here cj^ ^ and Cn,s are creation and annihilation operators for an electron of spin s on site n; m„ is the 
displacement of the n-th ion from its equilibrium position, so that — Un) is the deviation of the length 
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of the n-th bond from its equilibrium length. The first sum (lesc;ribes hopping with transfer depending 
linearly on bond length. The energy associated to cr-bond compressibility is described by the second term, 
K being an elastic spring constant. The hamiltonian of eq. (1) can be rescaled and rewritten in terms of 
the dimensionless coordinates /?„ = a(u„+i — Un)/to as 

T = ■ + cl,Cn+l,s) + 7 l^/3n (2) 

n,5 n 

where 7 = (Kto) / (2q^ ) accounts for the strength of the e-1 coupling (small 7 corresponds to strong coupling) . 
A diagonalization of the first term in the r.h.s. of eq. (2) gives the single particle electronic energy levels 
^m.siiPn})- For a given set of occupation numbers i'm,s, it is possible to determine the values of the 
coordinates /?„ {i.e., of the hopping constants) which minimize the total energy 

m,s i 

for a given 7. Here the first sum runs over the possible single particle energy levels; the second is a sum over 
the [N — 1) bonds. It is through such a minimization procedure that the models of ref. [2-3] account for e-J 
interactions. In the ground state of the half filled system this procedure leads to Peierls dimerization [10]. 

Within the framework outlined above, the first excited state is obtained by moving one of the two 
electrons occupying the A?'*'' level to the {N + 1)*'* level. In a long chain (even N 00) the set of bond 
lengths (i.e., the set of values of /3„) which minimizes ET{{Pn}) for this electronic configuration displays 
two kinks; these delimit a central portion of chain where the dimerization is inverted [10]. Corresponding to 
this bond geometry the N*^ and {N + 1)*^ levels are degenerate and are found at the center of the Peierls 
gap. Similar drastic distortions of the ground state bond geometry with inverted dimerization in the middle 
of the chain occur for short chains; however finite size effects modify the kink bond geometry and break 
the degeneracy of N*^ and {N + 1)*'* levels. In all cases the total energy corresponding to the optimized 
bond geometry is substantially smaller than the energy that the system would have for the same electronic 
configuration in the ground state bond geometry. From the Franck-Condon principle it should be expected 
that absorption experiments probe the situation where bond lengths are held to their ground state values 
while fluorescence experiments probe the spectrum found by optimizing the bond geometry. 

Similar considerations hold for higher electronic excited states. Indeed, for realistic values of 7, if 
one allows bond geometry relaxation, the lowest ^Ag excited state corresponds to moving both the iV*'' 
level electrons to the {N + 1)*^ level: in the N ^ 00 limit both the total energy and the bond geometry 
corresponding to this situation are the same as those of the first ^Bu excited state [12]. However, if the bond 
geometry of the ground state is kept fixed, the electronic configuration which gives the lowest excited ^Ag 
state, is diff'erent: it corresponds to moving one of the N*^ level electrons to the {N + 2)*^ level [10]. 

The rest of this letter is devoted to studying how adding Hubbard terms of the form 

r = + ^rn,t nm,l + Vi Tim Tlm+l (4) 

to 

m m 

to the hamiltonian of eq. (2) modifies the picture presented above. Here, as usual vo and vi describe on site 
and nearest neighbor e-e repulsion, is the number operator for electrons on site m and n^.T {nm,i) is the 
number of spin up (down) electrons. The distinguishing feature of the treatment presented here resides in the 
way the hamiltonian (sum of (2) and (4)) is dealt with. The fermionic part of the hamiltonian is diagonalized 
to give the (many body) energy levels Em{{Pn}, vo,vi). Then the set of values of the coordinates /3„ which 
minimize the total energy 

ET,m{{Pn}, Vo,Vi) = Em{{l3n}, Vo,Vi) + 7^13! (5) 
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associated to the; r/?"' (many body) level can be determined for given values of 7, fo and vi. This procedure 
is the natural extension of that of ref. [2-3]. The hopping constants (bond lengths) are not forced into 
configurations which cease to be optimal when e-e interactions are turned on. 

The treatment outlined above differs from the explanations usually given [13-15] to rationalize the 
spectroscopic results of ref. [4-9]: these are based on the results of Pariser -Parr-Pople quantum chemical 
calculations, where the hopping constant are forced into a dimerized configuration fixed from the outset. 
Within this scheme lattice relaxations are prevented: i.e., neither the hopping constants for the ground state 
nor those for the excited states are optimized. Bond lengths are obtained a posteriori from tt bond orders. 
Hayden and Mele [16] addressed the issue of geometry optimization in models including c c interactions: 
using an RG method they did obtain the optimized ground state geometry. However, they computed the 
energy of the excited states using the ground state geometry: this procedure does not account for the e-1 
effects underlying the soliton physics. 

The program described above has been implemented numerically within the full basis set of singlet 
states for half filled systems with up to eight carbon atoms [17]. The valence bond basis of ref. [18] for 
the S = subspace is used as a starting point. A symmetric fermion hamiltonian is obtained by changing 
to a new basis of singlet states by Gram-Schmidt orthogonalization; standard algorithms can then be used 
to find the required eigenvalues [19]. The set of coordinates which minimizes the r.h.s. of eq. (5) must 
correspond to bond geometries symmetric with respect to the midbond. Therefore for a system of N sites 
minimization of a function of N/2 independent variables is required: the downhill simplex method has been 
used for this purpose [20]. 

Figure 1 displays examples of results obtained in this way [21]: it shows the energy (relative to the 
ground state energy) of the 1 and of the 2 ^Ag states for two different values of 7 (7 = .9 [22] and 
7 = 1.2), N = 8 and vi = 0. The energies of l^B^ and ^^Ag obtained keeping the ground state bond 
geometry fixed are also shown. Note that level crossing [23] between 1 ^JB„ and 2 occurs at much lower 
values of vq for the optimized excited state bond geometries than for bond lengths fixed to their ground 
state values. Also larger 7's, e.g., smaller e-1 interactions, lead to l^Bu-2^Ag crossings at lower values of 
vq, at least for N < 8. Including a nearest neighbor interaction (non vanishing vi) does not (for reasonable 
values of the ratio v\/vq) change the qualitative features of these results. In the range < vq < 5, the 
main effect is to increase slightly the 2 ^Ag energy while that of 1 ^Bu is nearly unchanged. As a result, the 
l^Bu-2^Ag crossing occurs at slightly higher values of vq. Results qualitatively similar to these are found 
both for N = 6 and N = 4. A detailed description of these and of the other numerical results summarized 
in this letter will appear in a forthcoming publication. 

It should be noted that three dimensionless parameters (7, vq and vi) completely determine the ratios 
between the energies of the electronic states as well as the relative size of the hopping integrals. On the other 
hand in order to estimate bond lengths and absolute energies additional phenomenological constants [10] 
are needed. To avoid introducing other parameters, table 1 shows examples of how the hopping constants 
(rather than the bond lengths) change as e-e interactions are turned on: note that large hopping constants 
correspond to short bonds and vice versa. The experimental values of the energies for the 2^Ag and 1 
states are close in polyenes: i.e., realistic values of vq correspond to the 1 ^Bu - 2 ^Ag crossing region. It is 
clear from table 1 that the various types of soliton-like bond geometries (and in particular the reversed bond 
alternation in the chain center) survive at these levels of e-e repulsion. 

It seems appropriate at this point to comment on a recent letter by Konig and Stollhoff [24] which 
called into question the importance of the Peierls mechanism in determining the ground state dimerization 
of irans-polyacetylene. The results of table 1 for the half filled ground state are at variance with the 
conclusions reached by these authors. These results show that (for realistic [22] values of 7) e-e interactions 
have little effect on the ground state hopping constants and are consistent with a picture where the Peierls 
mechanism is the main reason for dimerization. The results of Konig and Stollhoff appear due to their failure 
to independently fit their ab initio results to those obtained from a semi-empirical hamiltonian which does 
not include correlations. 

In order to model systems where ground state degeneracy is lifted as a consequence of the local molecular 
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structure, an explicitely biased hopping term 



to 



to 



E (-1)" • 



) 



(6) 



has been added to the hamiltonian (here % is a site independent phenomenological constant). Figure 2 
shows numerical results for this situation when {tb/tg) — .08, 7 = .9 and vi — 0. As anticipated the 
1 ^Ag crossing occurs for higher values of vq than before. Again, analog behavior has been obtained 

for systems with = 4 and N = 6 and the distorted "bipolaron" -like bond geometry survives in presence 
of e-e repulsion. Although these results agree with the qualitative arguments presented at the beginning of 
this letter, in order to account quantitatively for the findings of ref. [6 7] computations on larger systems as 
well as more realistic forms of the terms lifting ground state degeneracy are needed. It should be stressed, 
in this regard, that the energies -Br.m computed in this letter refer to the semiclassical minima relative to 
th(^ optimized bond geometries for the m*'' many electrons level. Spectroscopic experiments, on the other 
hand, probe the various vibronic levels associated to this electronic state. 

The results of figs. 1 and 2 show that 1 ^Bu - 2 ^Ag crossing occurs for sufficiently high values of vq even 
if the bond geometries are held to their ground state configurations. However, failure to accomit for lattice 
relaxation for the electronic excited states [14] amounts to ignoring a physical ingredient which is essential 
in interpreting the available spectroscopic evidence. 

In summary, numerical results from a full many body description of conjugated chains which includes 
both e-J and e-e effects have been presented for chains with up to eight carbon sites. These systems are too 
small to allow reliable extrapolation to AT ^ 00 of detailed numerical results such as those for the energy 
of the excited states (for given 7, vo and vi) or for the strength of the e-e repulsion at which 1^3^-2^ Ag 
crossing occurs. However, in view of the results of ref. [10-11] it is natural to expect that for long chains the 
non linear excitations predicted on the basis of the models of ref. [2-3] will continue to correspond to the 
relaxed (minimum energy) bond geometries of the excited states when realistic e-e interactions are turned 
on. Moreover, the effect of such interactions on the electronic excitation spectrum for long chains [12] will 
be qualitatively similar to that discussed here for smaller systems. 

Several useful discussions with Ken Hass, Phil Pincus and Bill Schneider are gratefully acknowledged. 
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Table 1. 



Results for the hopping constants (1 — /3n) for chains with A'' = 8, 7 = .9 and vi = 0. The n*'* column 
gives the n*'' hopping constant with the outside bond corresponding to n = 1 and the central bond to n = 4. 
Ne is the number of electrons. The last set of data refers to a doped chain. 
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Figure captions. 

FIGURE 1. Energy Et (in units of to), relative to the ground state, of the 1 ^Bu state (empty squares 
7 = .9 , filled squares 7 = 1.2) and of the 2^Ag state (circles 7 = .9 , bullets 7 = 1.2). Here vi — and 
N = 8. The continuous (7 = 1.2) and broken (7 = .9) curves show the corresponding results obtained from 
the fixed ground state bond lengths. 

FIGURE 2. Same as in fig. 1 in a system where ground state degeneracy is broken by a biased hopping 
term of the form (6), with {th/to) = .08. Here 7 = .9, t;i = and N = 8. 



PostScript files containing figures are available on request from: rossi@scilab.srl.ford.com 
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